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We  began  research  on  developing  quantum  lattice  algorithms  in  the  Spring  of  2001. 
During  the  period  of  this  AFOSR  grant  we  have  tested  the  accuracy  of  our  quantum 
algorithms  on  exact  soliton  solutions  to  the  Nonlinear  Schrodinger  equation  (NLS)  as 
well  as  to  ID  MHD  turbulence.  The  following  papers  have  been  published  and/or 
submitted  for  publication  under  the  AFOSR  grant: 

Lattice  Boltzmann  and  Quantum  Lattice  Gas  Representations  of  One-Dimensional 

Magnetohydrodynamic  Turbulence 

L.  Vahala,  G.  Vahala  and  J.  Yepez 
Phys.  Lett.  A306,  227-234  (2003) 

Quantum  Lattice  Gas  Representation  of  Some  Classical  Solitons 
G.  Vahala,  J.  Yepez  and  L.  Vahala 
Phys.  Lett.  A310, 187-196  (2003) 


Quantum  lattice  gas  representation  for  vector  solitons 
G.  Vahala,  L.  Vahala,  and  J.  Yepez 
SPIE  Conf.  Proc.  5105, 273  -  281  (2003) 

Inelastic  Vector  Soliton  Collisions:  A  Quantum  Lattice  Gas  Representation 
G.  Vahala,  L.  Vahala  and  J.  Yepez 

Phil.  Trans..  Roy  Soc.  London  (accepted  for  publication,  2004) 

Quantum  lattice  representation  of  dark  solitons 
G.  Vahala,  L.  Vahala,  and  J.  Yepez 
SPIE  Conf.  Proc.  5436,  (to  be  published,  2004) 


As  insights  into  fluid  turbulence  has  been  gained  by  examining  simplified  nonlinear  one¬ 
dimensional  (ID)  models  as  Burgers  equation.  Since  MHD  turbulence  involves  more 
complex  structures,  Yanase  (1997)  and  others  have  developed  a  ID  magnetized 
generalization  of  the  Burgers  equation.  In  its  simplest  form,  this  self-consistent  model  is 
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where  the  fluid  velocity  a  =  and  the  magnetic  field  B=R(x,r)  z-  ju  is  the 

viscosity  and  t]  is  the  resistivity.  This  ID  MHD  model  reduces  to  the  regular  Burgers 
equation  in  the  limit  of  zero  magnetic  field  :  5^0.  However,  for  non-zero  magnetic 
fields  one  notes  that  there  is  a  self-consistent  magnetic  pressure  that  mitigates  to  reduce 
the  Burger  shock  formations  in  the  fluid  velocity.  Moreover  this  ID  model  exhibits 
rugged  invariants  as  the  full  MHD  equations  :  i.e.,  in  the  inviscid  limit  -♦0  the 

invariants  are  the  energy  and  cross  helicity  : 
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In  turbulence,  it  is  critical  to  accurately  compute  energy  transfer  rates, 
energy  transfer  from  the  magnetic  to  velocity  field  is  given  by 


with  the  kinetic  and  magnetic  energy  transfer  equations 
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For  ID  MHD,  the 
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To  develop  a  quantum  lattice  representation  for  this  ID  MHD  model  we  first  discretize 
the  spatial  dimension  and  introduce  2  qubits/scalar  field  at  each  node.  To  prepare  the 
initial  state  one  considers  each  qubit  to  have  the  form 

kaW)  =  4^  10  +  |0)  >  a  =  1...4 

where  is  the  probability  that  the  excited  state  is  occupied.  From  the  given  initial 
velocity  and  magnetic  field  profiles,  one  initializes  the  excited  state  probabilities 

Pi  =  i(l  +  M  +  B)  =  P3  ,  p2  =  ^(l+“-^)  =  Pa 

One  now  applies  a  local  unitary  collision  interaction  JJ  to  entangle  the  on-site  qubits 

|v'W)  -  , 

where  the  pie.collision  on-site  ket  |t/’(jc)J  =  |^i(a:  ))  ®  |?2(0)  ®  k3(0)  ®  k4(0)- 
On  these  post-collision  probabilities  one  performs  a  nonunitary  measurement  that 
destroys  the  quantum  entanglement  introduced  by  the  collision  operator: 

P'a{x)  = 

Finally,  one  streams  these  probabilities  p'„  to  nearest  neighbor  nodes. 

Thus  the  final  kinetic  equation  for  the  post-collision  probabilities 
p^{x+e^,t  +  l)  =  p„{x,t)+  {^p{x)\U^ h^U  -  n«k'(^)) 


To  recover  the  ID  MHD  equations,  one  first  entangles  qubits  ‘3’  and  ‘4’  with  the  square- 
root-of-swap  collision  operator 
/I  0  0  0\ 


U 


'swap 


0  Y 


l+i 

2 

0 


i+L 

2 

l-i 

2 

0 


and  then  streams  these  post-collision  qubit  states  ‘3’  and  ‘4’  to  nearby  sites 

+  s.k'W)  ■  ^.-(0.0.1, -I) 

One  now  reinitializes  the  occupation  probabilities 

-  pP  -  p^^' ,  pinp<‘>-pi"-pi'»  -  pf  -  pi^> 

and  then  entangles  qubits  ‘l’-‘2’  and  qubits  ‘3’-‘4’  followed  by  appropriate  streaming  to 
nearby  spatial  nodes.  Thus 


with  streaming  =  (l,-l,l,-l)  and  the  ;r/4  -  unitary  rotation  matrix 
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One  identifies  the  velocity  and  magnetic  fields  from  the  final  occupation  probabilities 
2u[x,t)  =  +  p^2\x,t)  +  p^i\x,t)  +  p^^\x,t) 

2B{x,t)  =  p^?{x,t)  +  p^}\x,t)  -  pf^{x,t)  -  p^i\x,t)  (4) 

and  the  unconditionally  stable  explicit  finite  difference  scheme  that  emerges  from  this 
quantum  lattice  algorithm  (in  lattice  units  Ax  =  1  =  At) 

«(x,t  +  l)  =  i[M(x-l,t)  +  M(x  +  l,t)]  +  ^^u\x-l,t)-u\x  +  l,t)^  + 

—[b{x -2,t)- B{x  +  2,t)]  [5(x - l,t)  +  2B(x,t)  +  B(x  +  l,t)] 
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B(x,t  +  l)  =  ^[B(x-2,t)  +  2B(x,t)  +  B(x  +  2,t)]  + 

— [m(x  -  l,t)B(x  -  2,t)+  {m(x  -  l,t)  -  u{x  +  l,t)}B(x,t)  -  u{x  +  l,t)B(x  +  2,t)] 


(5) 


In  the  continuum  limit,  Eq.  (5)  reduces  to  the  desired  ID  MHD  equations  (1)  with 

/a 

transport  coefficients  2a1  ’  This  finite  difference  scheme  can  be 

immediately  extended  from  a  scalar  to  vector  magnetic  field  B  -*•  B  =(0,Bj,,B^). 


A  typical  simulation  is  shown  below.  The  initial  oscillatory  velocity  field  (dashed 
curves)  steepen  following  a  Burgers-like  shock  development,  but  energy  is  transferred  to 
the  magnetic  field  so  that  these  field  are  amplified  in  the  regions  of  the  shocks  and  the 
velocity  shocks.  The  kinetic  energy  spectrum  exhibits  the  k~  spectrum. 


